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We study the linearized kinetic equation of relaxation model which was proposed by 
Bhatnagar, Gross and Krook'^'^^ (also called BGK model) and solve the dispersion relation. 
Using the solution of the dispersion relation, we analyze the relaxation of the macroscopic 
mode and kinetic mode. Since BGK model is not based on the expansion in the mean 
free path in contrast to the Chapman-Enskog expansion, the solution can describe accurate 
relaxation of initial disturbance with any wavelength. This non-relativistic analysis gives 
suggestions for our next work of relativistic analysis of relaxation. 

§1. INTRODUCTION 

Macroscopic phenomena of dissipative fluid are in general described by Navier- 
Stokes equations. However, Navier-Stokes equations can not provide good descrip- 
tions if the characteristic length of phenomenon is comparable to the mean free path 
of the constituent particles. For this reason Navier-Stokes equation cannot describe 
shock structure and short wavelength sound wave propagation^)"^-* that are impor- 
tant in astrophysics and plasma physics. 

To study short wavelength phenomena we have to use kinetic theory description, 
such as by the Boltzmann equation. Unfortunately it is quite difficult to solve the 
Boltzmann equation and various methods for obtaining approximate solutions have 
been developed. The most standard method is Chapman-Enskog expansion^) that 
was developed by Chapman, Enskog and Hilbert. This method begins with an 
expansion of distribution function in the mean free path of constituent particles and 
simplifies the results at a given order by introducing expressions from the previous 
order terms. This method derives the Euler equations at zeroth order and the Navier- 
Stokes equations at the first order and determines the transport coefficients that 
cannot be determined by Navier-Stokes equations. For this reason it is natural to 
expect that higher order approximations improve the description. The second and 
third order approximation was proposed by Burnett.^-*"'') However these equations 
provide little or no improvement in describing short wavelength phenomena. It is 
because the Chapman-Enskog expansion is a sort of iteration process and therefore 
inappropriate for a series that is not convergent. Modification of the Chapman- 
Enskog expansion has been proposed, for example, by Chen et a/.^^"^'') who avoid the 
recycling of lower order results. This extension of the procedures provided improved 
descriptions of short wavelength phenomena in the case of BGK model of kinetic 
theory, though this procedure cannot describe phenomena whose Knudsen number 
is close to or larger than unity. 

Another procedure is the moment method introduced by Grad.^^^'^^^ This pro- 
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cedure is a expansion of the distribution function not in the mean free path but us- 
ing Hermite polynomials and its 13-moment expansion reproduces the Navier-Stokes 
equation in the leading order. However this expansion is asymptotic rather than 
convergent so that the further development of expansion does not easily produce 
very accurate solutions in the short wavelength limit. 

In this paper, we derive the dispersion relation of BGK model^^^ of the lin- 
earized kinetic equation and analyze the relaxation process using the solutions of 
the dispersion relation. Since we study the kinetic model equation, the solution of 
the dispersion relation is correct in the domain of applicability of BGK model. The 
accuracy of BGK model has been studied in many papers: the experimental data for 
phase velocities and attenuation rate^^ indicates that BGK model is almost accurate 
in any Knudsen number^^)'^^-* but inadequate for phenomena sensitive to high-order 
moments. We are interested in time evolution of initial disturbance so that we 
solve the dispersion relation in terms of the frequency and obtain all modes of relax- 
ation. Since we use the solution of the dispersion relation of the kinetic equation, the 
solution accurately describes relaxation of initial disturbance with any wavelength 
and can be applied to various problem of the dilute gas dynamics; for example, appli- 
cations include dilute gas dynamics around an object,^^) study of the "atmospheres" 
of cosmic bodies, ^'^^ and the motion of a micrometer incident on the Earth's upper 
atmosphere.^^) In addition, this non-relativistic analysis may give suggestions to our 
next work of relativistic analysis of relaxation. 

In Section 2 we derive dispersion relations of BGK equation of the linearized 
kinetic equation. In Section 3 we solve the dispersion relation obtained in Section 
2 and derive eigenfunctions for all modes. In Section 4 we study the corresponding 
physical meaning of each eigenfunction. 

§2. THE LINEARIZED KINETIC EQUATIONS AND THE 
DISPERSION RELATION 



In this section, we derive dispersion relations of BGK model for the linearized 
kinetic equation. ^^^'^^^ The detailed calculation is presented in Appendix A. 
When there is no external field, BGK model of kinetic equation is 



Dfit, X, v) _ f{t, X, v) - /eg(i, X, v) 



(2-1) 



Dt T 

-=f# + v-vV (2-2) 
Dt \dt J ^ ' 

In Eq. (2-1), r is the relaxation time and /eg represents the local Maxwell-Boltzmann 
distribution 

P(t,x) (v-u(t,x))2 

/eg(t,X,v) = ,,3/2 ^"^^^^' (2-3) 

where 



R = kB/m, 



(2-4) 
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and /cs is the Boltzmann constant. 

Eq. (2-1) is a non-linear equation for /(t,x, v) because of the non-hnear depen- 
dence of feq on / through the following conditions: 



V'^ = m( 1, V, -(v-u)2 



(2-5) 
(2-6) 



These conditions are called matching conditions. 

For obtaining the dispersion relation, we start by expanding the distribution 
function around an equilibrium state /o(v) 

<5/ = /-/0, 5/e, = /e,-/0. (2-7) 

The linearized BGK model of the kinetic equation becomes the following form: 



Sf - if, 



We assume the following form of solution: 

Sf = Sfe~'''^^^~^°^~^^^'^. 



eq 



Then, Eq. (2-8) is 



where 



^-iuj + ik--v]Sf = ^5 feq, 



6 feq = fo 

/o 



^ + RTo ^ V2^ ~ 2j n 



Po 



(27ri?ro)3/2 



,-v^/2RTo 



m 



(2-8) 

(2-9) 
(2-10) 

(2-11) 
(2-12) 



Using the matching conditions, we can rewrite 6p, 6u, and 6T as the integrals 
of (5/ 



Su{t, x) 

(5r(t,x) 

Then Eq. (2-8) becomes 

- — ia; + • Sf(y) 
/o(v) 



6p{t,x) = J mSfd^v, 

m 



v5fcfv, 

Po 

2m ( 3, 



3po V2^ 2' 



To SfS-v. 



(2-13) 
(2-14) 

(2-15) 
(2-16) 



J ^ [po '^ToVWJh ' V2R^ ^ 3^ \2RTo 2 J [wJh vJ 
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In the following we take r as a unit of time, \/2RTq as a unit of velocity. 



UT = cj, \/2RTQTk = k, Z— 

Hereafter we omit". 

Finally, the linearized equation of BGK model is 

m{2RTo)^/^ /o(v) 



= V. 



(2-17) 



(2-18) 



where 



1 — ioj + ik-v 



1 + 2v ■ v' + 



./2 



(2-19) 



Above equations make sense only when 1 — iw + ik • v 7^ 0. We explain afterward 
the case where 1 — zo; + zk ■ v = 0. 

Eq. (2-18) is the homogeneous Predholm integral equation of the second kind. In 
particular, the kernel function K{v, v') can be separated with respect to the variables 
V and v' and this equation can be solved according to the general procedure. 

First we perform integration of Eq. (2-18) with respect to v. The equation 
becomes 

. ST ^ 
h37Tr = 0, 



(5/9 5u^ 



To 



where 



hi 



1 — iu 



V^e'''Erfc(6), 



/12 = 2i(l- 6Vvre^'Erfc(6) 



h3 

b 



b- (b^ + ^] V7Fe'''Erfc(6) 



1 — iuj 
k ■ 



(2-20) 

(2-21) 
(2-22) 

(2-23) 
(2-24) 



Secondly we multiply Eq. (2-18) by k • v and perform integration with respect 
to V. The equation becomes 



hi — + h2 
Po 



5Ur 



5T 



where 



721 = -i(l- 6V^e'''Erfc(6)) , 

722 = -k + 2b(l- 6V7re'''Erfc(6)) 



i23 



ib - (^6^ + i j 0Fe'''Erfc(6)| . 



(2-25) 

(2-26) 
(2-27) 

(2-28) 



Analysis of the relaxation process 



5 



Thirdly we multiply Eq. (2-18) by v± and perform integration with respect to 
V. The equation becomes 



I±±- 



6u I 



= 0, 



I±± = k- 0Fe^'Erfc(6). 



(2-29) 
(2-30) 



Finally we multiply Eq. (2-18) by v'^ and perform integration with respect to v. 
The equation becomes 



-131 r 132 

PO 



ST 



where 



I31 
I32 



^ [SA: - 2 {ft - {h^ - l)V^e'''Erfc(5)} 
- [26 {ft - (6^ - l)V^e'''Erfc(6)} - 3 

133 = \^^- {2^1 - b^) + C^b^ -b^ + l)\/^e'''Erfc(5)} 



(2-31) 

(2-32) 
(2-33) 
(2-34) 



If the determinant of the above homogeneous system is set equal to zero, the 
following dispersion relation is obtained: 



111 


/12 


/l3 








I21 


I22 


I23 








I31 


I32 


hs 

















I±± 




















= 0. 



(2-35) 



This condition implies either 



or 



111 


h2 


h3 


hi 


I22 


I23 


I31 


I32 


I33 




I±± 


= 0. 



0, 



Eq. (2-36) corresponds to longitudinal mode {5ux 7^ 0, 5u±_ 
corresponds to transverse mode {5ux = 0, 5u_\_ / 0). 

Finally we obtain the exact solution 5f in the following form: 



(2-36) 
(2-37) 

0), and Eq. (2-37) 



''^ ' ^ 1 - ?a;„ + zk • V 



+ 



2i?rn 



^1 .T. 



i(a;nt+k-x) 



(2-38) 



where C„ is a constant coefficient and (5u^j„, ^T^^ are eigenfunctions obtaining 

from Eq. (2-36). 
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We should study the case where 1 — + ik • v = 0. In this case the mode 
becomes continuous. According to Eq. (2-16) the eigenfunction for this mode satisfies 
the following equation: 







m m V 



2m ( 3\ / 



,,/2 



3/90 V2-R?o 2 ) 2 



(2-39) 



This mode represents decaying of the moments of / with vanishing (J/e^, i.e., = 
5T = 0,Su = 0. 

§3. RESULT 

3.1. DISPERSION RELA TION 

In this section wc show dispersion relations obtained in the previous section. 
First, we solve the dispersion relations numerically and results are shown in Fig. 1, 
Fig. 2, Fig. 3. 





Fig. 1. The decay rate of the thermal conduc- 
tion mode. 



Fig. 2. The decay rate of the shear flow mode. 





Fig. 3. The dispersion relation of the sound 
wave mode. The Solid line represents the 
decay rate (— Im w) and the dotted line rep- 
resents the frequency Re oj. 



Fig. 4. The dispersion relation of the trans- 
verse kinetic mode. The Solid line repre- 
sents the decay rate (— Im w) and the dot- 
ted line represents the frequency Re w. 
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Solid lines represent the decay rate (— Im u) and dotted lines represent the fre- 
quency Re CO. In the long wavelength part, solid lines are proportional to fc^ and 
reproduce the result that can be obtained by the first Chapman-Enskog approxi- 
mation. Regarding shear waves, we find only imaginary part of frequency and it 
indicates that shear waves cannot propagate in fluid when we consider initial value 
problems. 

100 I • • I 

3 

0.1 - 

0.01 r 

0.001 L^^j . — • — 

0.1 1 10 




Fig. 5. The dispersion relation of the longitu- 
dinal kinetic mode. The Solid line repre- 
sents the decay rate (— Im ui) and the dot- 
ted lino represents the frequency Re u). 



Fig. 6. Another dispersion relation of the lon- 
gitudinal kinetic mode. The Solid line rep- 
resents the decay rate (— Im uj) and the dot- 
ted line represents the frequency Re uj. 



We also find rapidly decaying modes — Im oj > 1/r in Fig. 5, Fig. 6, and Fig. 4. 
Since (— Im w) is 1/r at A; = 0, these modes represent kinetic modes that do not ap- 
pear in macroscopic descriptions. We discuss kinetic modes in the following section. 

Fig. 7 and Fig. 8 plot inverse of phase velocity and decay rates of sound wave 
modes that are obtained by solving Eq. (2-36) for A; as a function of given uj. These 
reproduce the results of Bhatnagar^^^ and are analogous to the results of Sirovich^^^ 
in which they have performed linear analysis of the Boltzmann equation for various 
forms of collision term. We show only sound wave mode for k but have found other 
modes, for example thermal wave modes. 
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Fig. 7. Phase velocity of sound wave mode in 
units of the adiabatic sound speed, ^J^RTo. 
The dispersion relation is compared to the 
experiment of Meyer and Sessler.^' 



Fig. 8. Decay rates of sound wave mode in 
units of the adiabatic sound speed, -^/^RTo- 
The dispersion relation is compaxed to the 
experiment of Meyer and Sessler.^^ 
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3.2. EIGENFUNCTIONS 

In this section we show eigenfunctions to study physical picture of each mode. 
Substituting dispersion relations obtained in Section 3.1 for Eq. (2-20), (2-25), (2-31), 
each mode is plotted with respect to the fluctuations (5p, 5ux and 5T that are nor- 
malized in 5p/po + 5ux/^(2RTq) + ST/Tq = 1. 




0.8 
0.6 
0.4 
0.2 






Sp 

6u 

ST . 







Fig. 9. Eigenfunction of thermal conduction 
mode. 



Fig. 10. Eigenfunction of sound wave mode. 
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Fig. 11. Eigenfunction of the longitudinal ki- 
netic mode. 



Fig. 12. Eigenfunction of the second longitu- 
dinal kinetic mode. 



In Fig. 9, ST is dominant and Su^ is very small for a long wavelength part, which 
indicates that this mode represents thermal conduction. Interestingly, Sux increases 
with decreasing wavelength and it indicates that picture of pure thermal conduction 
cannot be applied to the short wavelength part of this mode. 

In Fig. 10, Sp and 5ux are dominant for large wavelength and this mode can be 
regarded as sound wave mode. In the short wavelength regime ST dominates the 
others and it represents that this mode is oscillation between ST and Sux in this 
regime. 

Fig. 11 and Fig. 12 indicate that these two modes show similar relaxation. Both 
modes contain small Sp and large ST and Sux for any wavelengths. However, the 
timescale of decay of these modes is shorter than relaxation time r and these modes 
cannot propagate and do not appear on macroscopic scale. 
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§4. DISCUSSION 



4.1. KINETIC DECAY 



In general, relaxation is achieved by collisions so that decay rate (— Im uj) is not 
simply expected to be larger than 1/r. However, Fig. 3 and Fig. 2 indicate that it 
is contrary. In this section we discuss why this happens. 

In a short wavelength regime the relaxation rate (— Im cj) of this mode is larger 
than l/r and it contradicts our understanding that relaxation is achieved through 
collision. In the BGK model adopted in this paper, however, r is average relaxation 
time, namely r is independent of particle velocity v. Thus, we should consider that 
physical system become collisionless gas when k is larger. 

For obtaining physical description of relaxation in the collisionless regime, we 
consider macroscopic momentum equation. Prom kinetic theory, linearized momen- 
tum equation is 

Po^5u = -V5p-vV, (4-1) 

where P is traceless stress tensor. We consider only longitudinal mode. We multiply 
Eq (2-25) by oo and use ideal gas law: p = pRT. Then, Eq (2-25) can be cast into 
the following form: 

-k'^{6p + 6T)+u6p(-i + jio] (4-2) 



-ujk-du + k\dp + dT) + 

1 - 60Fe^'Erfc(6) + iojb - (^b'^ + ^ V^e^'Erfc(5) j | j 6T 



= 0. 



Fig. 13 plots the ratio of magnitude of each term of the momentum equation in 
sound wave mode Eq (4-2). This figure indicates that viscosity is small for a long 
wavelength and becomes larger for a short wavelength as in macro description. 




Fig. 13. The relative magnitude of each term Fig. 14. The relative magnitude of each term 
in the momentum equation for sound wave in the momentum equation for longitudinal 

mode. kinetic mode. 



Fig. 14 is the same plot for kinetic modes. These modes are purely kinetic and 
are not found in Chapman-Enskog approximation. Fig. 14 indicates that viscosity 
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dominates pressure and momentum density for a long wavelength. This is obvi- 
ously different from the description of viscosity of Chapman-Enskog approximation 
in which viscosity is small correction to macroscopic variables. This implies that 
these modes are purely kinetic modes. 

4.2. ASYMPTOTIC ANALYSIS OF DISPERSION RELATION 

For studying the property of the dispersion relation, we take the long wavelength 
limit.First we consider the shear flow mode. As in Sec. 2, the dispersion relation of 
the shear flow is 



m 



d\v±df 



m d vv± 



eg 



1 — iU! + ikVr 



5p V • (5u 

— H h 

Po RTo 



2RTa 



5T_ 



(4-3) 



We expand the integrand in powers of d = k/{l — iuj) and neglect the terms higher 
than second order on the right-hand side. Then above equation reduces to 

d'^v^^^ [1 - idv^ - d^vl] 6feg. (4-4) 



poSuj 



I 



Rewriting above equation, we get 

po5u^ = ^ ._ / d^yrv^{l - d^vl)5feq 

1 — liO J 

[1 — 2aid'^)poSu±. 



1 



1 — iu! 



where 



Neglecting u"^, we obtain 



ai 



2 2 

mvzv , 



Po{2RTo) 



2 feq- 



oj = —2iaid^ ~ —2iaik^. 



(4-5) 

(4-6) 

(4-7) 



Next we study the long wavelength limit of thermal conduction mode and sound 
wave mode. As in Sec. 2 the conservation of particle number is 



(5/9 = m 



Sf 



m d V- 



eq 



1 — ico + ikvx 



5p V • (5u 

— H h 

PO RTo 



2RTo 



3\ 5T 

2) n 



(4-8) 



The same as for the shear flow mode, we expand the integrand in powers of d = 
k/{l — iuj) and neglect the terms higher than second order on the right-hand side. 
Then above equation reduces to 



5p 



m 



1 — iu 
1 



j d\il - idVx - d^vl)5feq 

pod^ 6T 



dp — idpo 



V2mh 



2 To 



(4.9) 
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where we use 



I' 



(f\mvi5feq = pqRTq ( — + ^ ) . 



Rewriting the above equation, we get 

bp 



lu — 



id- 



= 0. 



2 7 po ^/2B% 2 To 
Similarly, the conservation of the momentum and energy reduce to 

id 5T 



id5p 2x ^'^x 

h (icj — zaAd ) , — — 

2 po V2RIh 2 To 



= 0, 



-iuj — 0:20?^ 



6p 

J 



2ia2d 



5u^ 



+ 



-iuj — d^ { as — -a2 



where 



02 



as 



04 



/ 



d\ 



2 2 
mv V 



— f 

n \9 Je 



po(2i?ro)2-'^'^ 4' 



f ,3,^^nvS^ _ 35 

J r 



d\ 



Pomnf" 8 
2 /eg 



(4-10) 
(4-11) 

(4-12) 
(4-13) 

(4-14) 
(4-15) 
(4-16) 



Po{2Rnf'"' 4 

Prom above equations we obtain the dispersion relation in the following forms: 



3^ . 
-2^ + 



19d^w 



2 ,-,2 



35 



+ I -id' + —id" \ ^- — 



5d^ 15^6 



16 



0. 



(4-17) 



We obtain three roots given by lot,cos± that are accurate to the second order in k 

id'^ iW 



W5± = Ml 



2 ' 
4z(i2 



(4-18) 
(4-19) 



We substitute d = k/{l — iu) into Eq. (4-19) and solve the equation with respect to 
u. Then we obtain 



(4-20) 



The above result shows that the decay rate (imaginary part) of the thermal con- 
duction Eq. (4-18) and the sound wave Eq. (4-20) is the same in the long wavelength 
limit. It corresponds to the fact that the value of the Prandtl number is unity in 
the BGK model. (The Prandtl number is the ratio of the viscosity times Cp to the 
thermal conductivity.) 
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4.3. MATHEMATICAL STRUCTURE 
4.3.1. DISCRETE SPECTRUM 

In Sec. 2 we have derived the hnearized equation of BGK model Eq. (2-18). 
However, as is indicated at the last part of Sec. 2, we assume 1 — + zk • v 7^ 
when we derived Eq. (2-18). 

When we assume 1 — itj+ik • v 7^ 0, we can formally think the linearized equation 
of BGK model, Eq. (2-18), as the following integral equation: 



'^/(v) = A J dVif(v,v')<5/A(vO, 



(4-21) 



where 



A = l, 
K(v,v') 



(4-22) 



m(2i?To)V2 /o(v) 



Po 



1 — zo; + zk • V 



1 + 2 V v' + 



,'2 



(4-23) 



This is the homogeneous Fredholm equation of the second kind with eigenvalue A = 1 
and obtaining frequency co is equivalent to deriving degenerated eigenfunction Sfx^j ■ 
Mathematically, discrete modes that include both macroscopic modes and kinetic 
modes are degenerated solutions of integral equation Eq. (4-21) and the number of 
discrete mode is equal to dimensions of eigenspace of eigenvalue A = 1. In Eq. (4-21), 
K{v, v') is the compact integral operator. This indicates that the number of discrete 
modes is not infinite but finite because of the Fredholm alternative theorem. 

4.3.2. CONTINUOUS SPECTRUM 

In Sec. 2 we obtained continuous spectrum 1 — ■io; + zk • v = in addition to the 
discrete modes. In this section we discuss the eigenfunction of this mode. 

We assume a situation that a beam of particles of constant velocity Vq is imposed 
in the initial condition. Then Sf becomes as follows: 



Sf = 6h + DSD{y-yo), 



(4-24) 



where D is constant coefficient and Sd is the Dirac delta function. D5d{^ — vq) rep- 
resents particle flux distribution of v = Vq. represents the perturbed distribution 
function of gas. 

We substitute this 5f into Eq. (2-16) then it becomes 



- — zoj + ik • V j 5fi + f - — za; + ik • Vq ) -D 



I /o(v 



3,,//o(v) 



m m V 
— + 



Po-RTo ■ 3po \2RTq 
m m V 2m f 3\ / 

Jo + JoRTo ■ ^ 3^ V2^ ~2) \2RTq 



2 J \2RTo 



_ 3 
~ 2 

D. 



(4-25) 
<^/i(v') 
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Then we assume next continuous spectrum. 



oj = h k • Vq. 

r 

We substitute this oj into Eq. (4-25). Then Eq. (4-25) becomes 



Sfi 



I 



v7o(v) 



m m V 
— + 



+ 



irk • (v - vo) [po po RTq 

/o(v) 



, 2m 
• v' + — 



3po \2RTo 2 J \2RT^ 



,,'2 



(4-26) 

(4-27) 



irk • (v — Vo) 



m ^ m V 
Po Po RTo 



Vo + ^ 







( vl 




\2RTo 


2T J 


\2RTo 





D. 



The solution of this integral equation can be obtained the same as in Sec. 2. 
We explain the case of shear flow for simplicity. We multiply Eq. (4-27) by v_\_ and 
perform integration with v. The equation becomes 



I A; - V7re'''Erfc(6)} Su^y = V7re^'Erfc(6) D voy, 

ik • Vo 



k 



-IVOx- 



Finally we obtain eigenfunction 5uiy in the following forms: 

-1 



5u 



12/ 



k 



1 



DVOy. 



(4-28) 
(4-29) 

(4-30) 



5pi, (5ui, 5Ti can be obtained with similar procedure. 

4.3.3. GENERAL SOLUTION 

The continuous mode and the discrete mode we discussed in previous sections 
seem to cover all spectrum of Eq. (4-21) and the general solution of the linearized 
BGK equation would be as follows: 

<^/(v) = ^C„i^(v,v')(^/a,„(v') + / d3voC(vo)K(v, v')(5/.o(vO, (4-31) 

k{v, v') = {l-iu + ik- v)K{v, v'). (4-32) 

Because of the Fredholm alternative theorem the number of the discrete modes 
is finite. We may conclude that most of the degrees of freedom of 6f belong to the 
continuous spectrum and only finite parts of them belong to the discrete modes. 



§5. CONCLUSION 



We have solved linearized kinetic equation of BGK model and obtain the exact 
solution for the relaxation of initial disturbances that describe the kinetic modes in 
addition to macroscopic modes that can be obtained by Chapman-Enskog expan- 
sion. Using those solutions, we have obtained the eigenfunction and analyzed the 
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relaxation process corresponding to those solutions. Since we use the solution of the 
dispersion relation of the kinetic equation, the solution describes accurate relaxation 
of initial disturbance with any wavelength. In next paper we will apply this method 
to the relativistic kinetic equation. 
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Appendix A 

THE DETAILED CALCULATION OF THE FLUID EQUATIONS 



In Sec. 2 we obtained BGK model of linearised kinetic equation in the following 
forms: 

<5/(v)= j d\'K{-vy)6f{v'), (A-1) 



i^(v,v') = 



m(2i?To)3/2 /o(v) ^ 2vv' + ^ ^ ^ -'^ ^ 



Pq 1 — ioj + ik-v 



, V ] \ V 

3 V 2M 2 



(A-2) 



where /o is the unperturbed state. 

We study with the spherical coordinates in velocity space setting k as the z- 
axis. Then we have to relate v • c^u in Sfeq to k. We define 0, 6', & as the angle 
measuring between k and v, k and Su, v and Su, and ^ as the angle between v and 
6u measuring on the plane projected perpendicular to k. Then we have the following 
relation 

cos O = cos cos 6' + sin 6 sin 9' cos (/). (A-3) 

In this section we calculate the momentum equation of k direction as an example. 
We set X and y axis in the direction of k and Su± that is the projection of Su on the 
plane perpendicular to k. Then Sux becomes 

kSux = k • / d^vmvSf 



= 27rm Hy'^dv C d^ ''^'^^"^ , (A-4) 

Jo J-i + ikvn 

where jx = cos 9. Then we replace 5feq with Eq. (2-11) and replace dimensional 
variables with dimensionless form as in Eq. (2-17) 



dv r da I (v^ -l]6T + 5p + 2Tk • <^u ^ (A-S) 



t-v/tt Jo J-1 1 — ico + ikvp, \ \ 2 ) k 
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We remove n of odd index and the equation becomes 



Vt^t Jo 

4/90 



Jo 

where 



dv I dfi——-p, 7^ 

dv / d/j, 
Jo 



-ivfi \ v~ — - j ST — ivfxSp + 26— k • 



V e 



-I V 



•2 -]5T -i5p+^\i-5vi\ , 



2b, 



1 — iu; 
k ' 



Then we use next integral formula 



1 1 b 

I "A* ^5 PT^ = -~r arctan -, 



(A-6) 

(A-7) 

(A-8) 



and the equation can be integrated 



4po 
4po 



V^r Jo 



dvv^e-"' 



dve-" 



1 - - tan"^ -] \ -i{v'^ --]6T-i6p + —k-Su 



26. 



2 _ . 26, 



-i ( f ^ - - 1 (5r - z(5p + — k ■ du 



— 6tan ^ ^ 1^ (^~^^ + v — V + dTivSp + • Su 
Using the identity 



_d_ 
dv 



(A-9) 



(A-10) 



we perform partial integration 



4po 
-v/ttt 



/"OO /*! 

7o Jo 



62 + z;2 



•2 (5T-i(^p + ^k-(^u 



t;2 1 



JT + -Sp -^k-Su 
2 A; 



We use next integral formulas 



f 

Jo 



x2 + 62 26 
(o > 0, 6 > 0) 



Erfc(a6), 



J' 

Jo 



(a > 0) 



(2n-l)!! rH 



2n+l y Q2n+1 



(A-11) 

(A-12) 

(A- 13) 
(A-M) 
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Finally the equation can be integrated as follows: 



Po 

T 



ib 



6 - + V7Fe'''Erfc(6)| ST -i(^Sp + • Juj (l - 6 V^Fe''' Erfc(6)) 

(A-15) 

On the other hand, integration with respect to 6feq is clearly 

k • J dvmv5feq = ^k • Su. (A-16) 
In the end, the momentum equation is 



-k • (5u + ib 



ST 



-i (^Sp + ^k ■ Su^ (l - 6V^e'''Erfc(6)) = 0. (A-17) 

The energy equation can be performed in the same way. 

The shear flow equation can be calculated as follows. Su± is given by 

poSu± = J d^vmv±5f 

= m[d\v^ — Iv'^ --]ST + 5p + 2^-Su 

J 1 — iuj + ikvx L \ 2/ 

= m I d^Y •^"'^ , 2vlSu^. (A-18) 

J 1 — iu) + ikvx 

Then after the same calculation of the momentum equation we obtain Eq. (2-37). 

Appendix B 

ANALYTIC CONTINUATION 



Eq. (A-12) can be used only when 6 > and we have to regard b as if b has 

imaginary part, wc have to replace b to —b in Eq. (A-12). However in this paper 
we perform analytic continuation and remove these restrictions. In this section we 
explain the analytic continuation. 

For simplicity we explain by using the transverse shear flow mode {Su± ^0). In 
Eq. (A-18) we consider the circular cylindrical coordinates {v±, (f>, Vx) and perform 
(f) integral. Replacing to dimensionless form, the equation reduces to 



/OO P 
dvx 
-OO Jo 



'i{kvx — i{oj + i)} 
Performing v±_ integral the equation becomes 



OO „.3 

dv±-r-. ^, -ye-'''-^'. (B-1) 



1 l"^ e~'"^ 

7^ / dvx^, ^. (B-2) 
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This integrand may have first-order pole on the real axis. However this singu- 
larity is removed to the upper half-plane if Im w > —1 and this results from the 
collision term of BGK model Eq. (2-1). This fulfills the demand of the Landau 
method that this problem has to be considered as the Cauchy problem. For this 
reason the integration path has to be distorted below the singularity. 

Concerning above equation the restriction Re 6 > is not satisfied when Im w < 
— 1, so analytic continuation has to be performed. The equation has first-order pole 
so we use the principal integral 

The integration without analytic continuation corresponds to the integration 
path distorted above the singularity when Re 6 < 0, so the principal integral is as 
follows: 



e(-'')^Erfc(-6) = ip H dv^^—^ ^ - ^e-(^'^\ {B-A) 

^ ' 2 i{kvo, - i{oj + i)) 2ik ^ ' 



2k i-b) 

In the end, Eq. (B-2) is as follows when Re 6 < 



oo 

dvx- 

oo 



i{kvx — i{uj + i)} 



A4e(-^)^Erfc(-6) + ^e-(^»A + ^e-(*) 
2k -b ^ ^ 2ik J 2ik 



2 



= ^e^^(2-Erfc(-6)) 

= ^e^^Erfclfe). (B.5) 

This indicates that analytic continuation permits us to use Eq. (2-37) even when 
Re 6 < 0. 
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